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fN| Abstract 

j>! The renormalized entanglement entropy (REE) across a circle of radius R has been pro- 

posed as a c-function in Poincare invariant (2 + l)-dimensional field theory. A proof has been 
presented of its monotonic behavior as a function of R, based on the strong subadditivity of 
entanglement entropy. However, this proof does not directly establish stationarity of REE 
at conformal fixed points of the renormalization group. In this note we study the REE for 
the free massive scalar field theory near the UV fixed point described by a massless scalar. 
Our numerical calculation indicates that the REE is not stationary at the UV fixed point. 
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1 Introduction 

A well-known problem in quantum field theory is how to define a quantity which decreases 
along any renormalization group (RG) trajectory. In two-dimensional QFT, an elegant 
solution to this problem was given by Zamolodchikov IT], who used the two-point functions 
of the stress-energy tensor to define a monotonic "c-function." At RG fixed points the 
Zamolodchikov c-function equals the Weyl anomaly coefficient c. 

An important additional property of the Zamolodchikov c-function is that it is stationary 
at the fixed points. For example, if we consider perturbing a CFT by a slightly relevant 
operator O of dimension 2 — 5 then 

c(g) = c uv - g 2 6 + 0(g 3 ) , (1.1) 

where g is the renormalized dimensionless coupling. More generally, in a theory with more 
than one coupling constant, 

% = , (1-2) 

where Gy is the Zamolodchikov metric, and /3 J = jj,^- is the beta-function for the coupling gi . 
The fact that the metric is non-singular guarantees the stationarity of the Zamolodchikov 
c-function at any fixed point in two-dimensions. 

It is of great interest to find out if these results extend to field theory in dimension d > 2. 
In four-dimensional conformal field theory there are two Weyl anomaly coefficients, a and 
c. Long ago Cardy conjectured [2] that it should be the a-coefficient that decreases under 
RG flow. This coefficient can be calculated from the expectation of value of the trace of the 
stress-energy tensor in the Euclidean theory on S 4 . Using conformal perturbation theory, it 
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is possible to establish the analogue of (1.1) in four-dimensions (2] 



c(g) = auv - 9 2 S + 0{g 3 ) , 



;i.3) 



where the perturbing operator O has dimension 4 — 5. Over the years considerable evi- 
dence was accumulated in favor of Cardy's conjecture, especially in supersymmetric 4-d field 
theories where a is determined by the U{1)r charges [3]. In particular, the principle of a- 
maximization |2), which states that at superconformal fixed points the correct R-symmetry 
locally maximizes a, has passed many consistency checks that rely both on the field theo- 
retic methods and on the AdS/CFT correspondence [B]-j7j. Last year, a general proof of the 
a-theorem was constructed in |8j[9]. 

A proposal for a similar theorem, called the F-theorem, in 3-d Euclidean CFT was made 



in 10 : 



log \Z S z\ 



;i-4) 



where Z$z is the Euclidean path integral of the CFT conformally mapped to S 3 . The 3-d 
analogue of a-maximization is that the R-symmetry of M = 2 superconformal theories in 



three-dimensions locally maximizes F 11,12 . In 10 it was further conjectured that in 



unitary 3-d CFT F is positive, that it decreases along any RG flow, and that it is stationary 
at the fixed points. The three-sphere free energy is ambiguous along RG flow but is well- 
defined for any conformal field theory. Thus, a precise statement of the conjecture is that, if 
there is smooth RG flow from a UV CFT to an IR CFT, then their three-sphere free energies 
satisfy Fuv > -Fir- Various Af = 2 supersymmetric examples have been presented in support 



of these statements; see, for example, 10 , 13 , 14 . Perturbations of a CFT on a three-sphere 



by slightly relevant operators were studied in 15 , and results analogous to ( 1.1 ) were found 



In particular, the S 3 free energy (1.4) is stationary at the fixed points. 



on 



Closely related ideas in Minkowski signature were advanced in 16 . For a field theory 
IR 2 ' 1 , it was proposed in fl6) that one should consider the finite part of the entanglement 



entropy S(R) across a circle of radius R. It was further shown in [XT] (see also [18]) that if 
the field theory in question is a CFT, this entanglement entropy agrees with the free energy 
of the Euclidean CFT on S 3 : if such a field theory is conformal, then 



S(R) 



2-kR 



a- 



1.5) 
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where e is the short-distance cutoff. 

The ideas related to the quantum entanglement entropy were recently made more precise, 
and a proposal for a monotonic c-function for a general, Poincare invariant (2+l)-dimensional 



QFT has emerged. It is the so-called Renormalized Entanglement Entropy (REE) [19 20 



f(R) = S(R) + RS'(R) . (1.6) 

An advantage of this definition is that the UV divergence proportional to the radius cancels, 
so that REE is manifestly finite for theories that are conformal in the UV. Studies of various 



explicit examples, both in free field theory and in the AdS/CFT correspondence 19 21 
yielded a monotonic c-function JF. Furthermore, J-'(R) = RS"(R). It was shown in [20] 
using Strong Subadditivity of the EE, that for any Poincare invariant field theory S"(R) < 0. 
This demonstrates that F{R) is a non-increasing function and establishes the F-theorem. 

In this note we address a question about the proposed c-function J-'(R) that has not 
yet been elucidated: namely, is it stationary for arbitrary perturbations around a CFT? In 



studies of holographic entanglement entropy 22,23 , T turned out to be stationary 19,21 



but it is not clear that this is a general property of the REE. 

In a general field theory which does not have a gravity dual, the calculation of REE 
is a difficult problem even if we resort to numerical methods. In this paper we consider a 
particularly simple example of RG flow provided by a free massive scalar field on IR 2,1 with 
the action 

I = -\f A[(W 2 + m¥] . (1.7) 



In this case there are efficient numerical algorithms for calculating T 24 ,25 . On dimensional 
grounds, T must be a function of g = (mR) 2 , which is the dimensionless coupling associated 
with the relevant perturbation m 2 (f) 2 /2. The UV fixed point of this theory is that of a 



massless free scalar, for which the value of J 7 is known analytically 15 , 18 ; therefore, 



^0/ = 0)=Jbv = ^- ^« 0.0638 . (1.8) 

8 167T^ 

The stationarity of the c-function at the UV fixed point would require ^-"'(O) = 0. However, 
our numerical calculation instead indicates that F'ig) is negative at small g. In fact, it may 
diverge in the limit g — > 0, but the limitations of our numerical work do not allow us to 
make a precise statement about the small g behavior of F'ig)- 

Our numerical results suggest that REE does not in general define a c-function in the 
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Zamolodchikov sense, because it is not always stationary at conformal fixed pointsj^] It is 
not clear how to compare this statement with the proof 20 of the monotonicity of J 7 . If g 
can have either sign, then ^'(0) 7^ would violate monotonicity. In the example we have 
considered, however, insisting on a stable UV fixed point requires that g is positive. So, 



there is no immediate conflict with the work of 26 . We note, however, that the absence of 



stationarity means that an equation like (1.2) cannot in general apply to T. Examination 
of other examples where the stationarity of T can be tested is highly desirable. 



2 Strip entanglement entropy at small mass 

Although the renormalized entanglement entropy across a circle is of our main interest, as 
a warm-up we start with a slightly different but simpler quantity, the entanglement entropy 
of a strip. 

The entanglement entropy between the strip of width R and its complement in (2 + 1)- 
dimensional free field theory is simply related to the entanglement entropy between the 
interval of length R and its complement in the (1 + l)-dimensional theory. More precisely, 



for the (2 + l)-dimensional massive scalar field of mass m we have the relation 27 



sffV.*) = \ rdpSS^(y/F+^,R) , (2.1) 
» Jo 

where \/p 2 + m 2 is the mass of the (1 + l)-dimensional scalar field, and L is the length of 
the strip. To derive this equation one first compactifies the direction parallel to the strip 
to a large circle of length L. Decomposing the (2 + l)-dimensional scalar field into angular 



momentum modes along this circle and then taking L — > 00 leads to (2.1). 
It is useful to define the entropic c-function for the strip through 



CUrip = R? d R Sj£f^ , (2.2) 
where = ^iSp /-^ is the entanglement entropy per unit length. This function is 



manifestly finite and cut-off independent 28 , 29 . For the massive scalar field, C str i P is also 



simply related to the (1 + l)-dimensional entropic c-function, defined 26 by c(t = mR) 



1 In (1 + l)-dimensions the monotonic c-function derived from the entanglement entropy of a segment |26| 
is not stationary with respect to m 2 either. In (1 + l)-dimcnsions this is likely due to the fact that cf> 2 is not 
a conformal primary field in the CFT of a massless scalar. 
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R9rS££L, through 02.1): 



C s t r ip (mi?) 



2tt 



dx c(a/ x 2 + (mR) 2 ) 



(2.3) 



We should note, however, that in general C strip is not expected to be a good c-function; for 
example, it is not constant along lines of fixed points. 

The entropic c-function for the (1 + l)-dimensional massive scalar field is a well studied 
quantity (see [30] for a review). The function interpolates between 1/3 in the UV (t = 0) and 
in the IR (t — > oo). At small and large values of t one can calculate the leading behavior 



of c{t) analytically 27 



c(t) 



1 1 
3 + 21ogt + 

c(t)= t l K 1 {2t) + ..- 



t < 1 
t > 1. 



(2.4) 



Unfortunately there is no known closed form expression for c(t) along the entire RG flow. 
The function may be constructed numerically by solving an infinite sequence of nonlinear 
differential equations 27 . Following this procedure leads to the curve in figure [TJ 




Figure 1: The entropic c-function c = R8r S^^ al for the (1 + l)-dimensional massive scalar 
field as a function of t = mR, where m is the mass and R is the length of the interval. 



The black curve comes from a numerical calculation using the prescription in 27 . The blue 



and orange curves are the analytic approximations in (2.4) at small and large values of t, 
respectively. The dotted red line marks the conformal value c(0) = 1/3. 



The function C s tn P (mR) is a monotonically decreasing function that approaches zero 
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exponentially fast in the IR. Using (2.3) one may determine numerically its value at mR 



27 



Cstrip(O) 



1 f°° 

- / dtcCt) w 3.97 x 10" 
* Jo 



(2.5) 



At large mR we may use the approximation for c(t) in the second line of (2.4) to write 



C s t r ip (mR) « — e 



mi? » 1 



(2.6) 



In figure [2] we numerically plot C st rip along the RG flow together with the IR analytic ap- 
proximation. 



0.04 



strip 




Figure 2: The function C str i P = R 2 8r S^l for the free massive scalar field in black, where 
^strip^ = ^stripV-^ * s the entanglement entropy per unit length across the strip of width R. 



The orange curve is the IR approximation in (2.6). The dotted red line is the initial value 



at t — mR = given in (2.5). 



It is interesting to ask whether C s t r i P is stationary at mR = 0. The first derivative of the 



entropic c-function (2.3) with respect to t — mR gives 

t 



^strip(^) 



1 

27 
1 

2^ 

->• 



dx 



— oo 
oo 



Vx 2 + t 2 
1 



z'(Vx 2 + t 2 ) 



-oo Vx 2 + 1 
1 



1 



2\Ac 2 + 1 log 2 {Wx 2 + r 



+ 0(log- 3 t) 



(2.7) 



t-+o 4 log 2 1 ' 
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where we used the asymptotic form of the (1 + l)-dimensional entropic c-function (2.4) 
and rescaled x — > tx in the second line. While this implies that Cg trip (0) = 0, stationarity 
additionally requires that C ; stri (t)/t vanishes at t — 0. This is clearly not the case since this 
quantity diverges as -l/(4tlog 2 t) as t -»■ 0. As can be seen in figure |3j we confirm this 
behavior numerically by plotting C' stTip (t)/t along with the analytic prediction. 




Figure 3: The first derivative C' stTi (t)/t as t = mR —>• 0. The black curve comes from 
the numerical calculation, and the dotted orange curve comes from fitting the numerics to 
the function C' stlip (t)/t ~ a/(t\og 2 t) as t — > 0. We find a ~ —0.25, which agrees with the 
analytic result in ( 2.7[ ). This means that C str i P is not stationary for the massive scalar field. 



To summarize the results of this section, we have found that the small mass expansion 



of the UV finite function (2.2) for a strip of width R is 



C 8tlip (mR) « 0.0397 - - 1 - + • • • . (2.8) 

4 log {mR) 

In the next section we will present evidence that the REE across a circle of radius R has a 
similar structure, with dJ^/dg looking clearly negative at small g = (mR) 2 . 



3 Disk entanglement entropy at small mass 

The entanglement entropy for the massive scalar field across the circle of radius R may be 
calculated numerically following the prescription in [24| , [25] , [30|[3T] . This numerical method 
has passed many non-trivial checks. For example, in [25] the terms proportional to (mR) 
and l/(mR) in the large (mR) expansion of the disk entanglement entropy were matched 
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numerically to the analytic predictions to high accuracy. A similar analysis in 32 matched 
the l/(mR) 3 term in the large mass expansion of the EE to the analytic prediction, and 



in [19] the value of the REE at mR = was shown to agree with (1.8). 



The numerical procedure 24 , 25 30 , 31 expands the scalar field into modes of integer 



angular momentum n. The radial direction is discretized into a lattice of N units. For each 
n the discrete Hamiltonian takes the form H„ - 



2 



K%(j)j, where 7Tj is the 



conjugate momentum to fa and i 
are 



N. The nonzero entries of the N x N matrix K n 



3 
2 



K? = "+n 2 + m 2 



71 2 

2 + — + m 2 



K 



Li+l 



K 



i + 1/2 



(3-1) 



To compute the entanglement entropy we need to know the two-point correlators X = 
(fafa) = \(K~ l l 2 )ij and P = {ni^j) = \{K +1 / 2 )ij. If the radius of the entangling circle R is 
a half-integer in units of the lattice spacing, then we must reduce the matrices Xij and Pij 
to the r x r matrices and P[j, which are defined by taking 1 < i,j < r with r — R — | . 
The entanglement entropy is then given by 



with 



Sn 



tr 



S(R) = S + 2j2S n 



71=1 



V^i^ + log (v 7 ^ + \ 



X rp r _ 

n n 



log (^/XrPr 



(3.2) 



(3.3) 



Further details of the numerical calculation can be found in 25 



In order to achieve the continuum limit we need to take N ^> r ^> 1 and m < 1 (restoring 
the lattice spacing, the latter condition is me 1). Then mR can be either small or large, 
and we can explore the REE as a function of this parameter. 

In our calculations we use a radial lattice consisting of iV = 200 points. We want to 
calculate the entanglement entropy for .06 < mR < 2, but to minimize lattice effects we 
restrict 30 < r < 50 in lattice units. To accomplish this we calculate the entanglement 
entropy for m = .002 • i in inverse lattice units, with i — 1, ... ,20. To take into account 



finite lattice effects, we follow 19 and for the lowest 10 angular momentum modes we repeat 
the calculation on lattices of sizes iV = 200 + 10 ■ j, with j — 0, . . . , 50, and then extrapolate 
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to iV = oo. We perform the numerical calculation for the first 3000 angular momentum 
modes. We take into account higher angular momentum modes by using the asymptotic 
behavior of S n derived in Appendix \A\ 

From the entanglement entropy we construct the renormalized entanglement entropy J 7 
along the RG flow. A plot of this function versus (mR) 2 is given in figure [4l From this plot we 



r 




Figure 4: The renormalized entanglement entropy T across the circle of radius R for the 
massive real free scalar plotted in black as a function of (mR) 2 . In this plot it can clearly 
be seen that <9( m ^)2 J 7 is negative and nonzero at (mR) 2 = 0, which implies that the REE J 7 
is not stationary at the UV fixed point of a free massless scalar field. The dotted red line is 
the zero mass value J-yv = ^ — ll^ - 

can see that d^^J 7 is clearly negative and nonzero at (mR) 2 = 0. Moreover, it appears that 
this quantity might actually diverge as (mR) 2 — > 0, though our limited numerical precision 
prevents us from making any precise statements along these lines. 



4 A toy model: the massive scalar on H 2 x S* 1 



For a CFT in (2 + l)-dimensions the calculation of the Renyi g-entropy of a disk may be 



mapped to a calculation on H 2 x 5 1 17,33 , where q is the radius of the circle. For a free 



scalar field of mass m the free energy on this space is [34 



7 q (m 



dXD(X) [log ( 1 - e ^W x +( mR A + nq^X + (mR) 2 



(4.i; 
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where the density of states is given by 



V(X)d\ = YglE) tanh(7rv / A)rfA 

4:71 

and the regularized volume of the hyperbolic space is Vol(H 2 ) 
finite part of entanglement entropy is 



(4.2) 



-27r. The formula for the 



Si 



dq 



J=l 



(4.3) 



q=l 



Applying this relation at m = 0, which corresponds to a massless scalar on a disk, one 
obtains Si = — + ?S^t, in agreement with other methods 



34 



The status of the calculation on M 2 x S 1 for m 2 > is less clear since this does not 
correspond to turning on mass on the original disk geometry. It is still interesting to inquire 
whether Si defined above is stationary with respect to turning on m 2 . An explicit calculation 
yields 



dSi 



dm 2 



<m?=0 



7T 

16 J 



R 2 



(4.4) 



Similarly, the calculation of the Renyi entropy can be mapped to that on the g-fold covering 



of 5* 3 34 . For a massive free scalar field of mass m on this space, the derivative of Si with 
respect to the mass of the scalar field is T^-| m 2 =0 = — f|-R 2 , which agrees in absolute value 



with (4.4) but has a different sign. Either way, Si is not stationary in these toy models. 
This lack of stationarity of Si, which is easily established analytically, is reminiscent of the 
numerical result for the disk entanglement found in the previous section. 
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A Large n asymptotics of S n 

A large n expansion of S n is made possible by the fact that to leading order in n the 



matrix K% defined in (3.1) is diagonal. Therefore, to leading order in n, the matrices 
X n = \K~ l/2 and P n = \K l J 2 are also diagonal, and so are their truncated counterparts X r n 



and Consequently, in this approximation X^P^ = j, and from (3.3) we conclude that 



S n vanishes. The non-trivial dependence of S n on n comes from the 1/n corrections to the 
relation X'^P^ = 1/4, which we will now calculate. 
From P 2 = jK n , we find for i > 2 



n 

pi,i+l _ pi+l 



n i{m 2 + 2) z 3 (m 4 + 4m 2 + 6) ^, 

h 7, + 0(n~° • 

2i An " 



lQn 3 
+ 11 ^ i 3 / 2 (i 



/ ' -(% + l) 3 / 2 (2 + m 2 



P 



i+2,i 



An 

+ 2) 3 / 2 



+ 0{n 



16n 3 



+ 0{n 



(A.l) 



In general, P^ l+a = 0(l/n 2a ~ l ), so up to a fixed order in 1/n the matrix P n is almost 
diagonal. From the relation X n P n — |, we can obtain a similar expansion for the entries of 
X„: 



\rii 



i i 3 (m 2 + 2) i 3 (A + 3i 2 (m 4 + Am 2 + Q)) 7 

n + — ^77^ + °( n ~ 

2n " 



X, 



An 3 16n 5 
i 3 ' 2 (i + l) 3 ' 2 % 3 I\% + l) 3 / 2 (3z 2 + 3i + 1) (m 2 + 2) 



An 3 



+ 0(n 



-7\ 



(A.2) 



i>n 



x M +2 _ x ^ _ ^ 2 (^ + 2p (3.2 + 6^ + 2) 

X n -X n — + U[n ) 



with the general asymptotic form X^ 1+a = 0(l/n 2a+1 ). By definition, X n P n = 1/4, or 
explicitly 



(X n P n Y' l+a = J2^n i+b P, 



■i+b,i+a _ ^aO 
n 4 



(A.3) 



The product of the reduced matrices X r n and P r , however, fails to equal 1/4 precisely because 
some of the terms in the sum over b are absent, namely the ones with i + b > r. Using 



the asymptotic expansions (A.1)-(A.2), we find that up to 0(n ) corrections, the matrix 
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X r n P^ — I is identically zero except for its lower-right two-by-two block: 



_ 1 r 2 (r + l) 2 _ (m 2 + 2)r 2 (r + l) 2 (2r + l) 2 
1 n nJ ~4 + 16n 4 " 32n 6 +t7ln > 



, n,; r.r-1 ^(r - lfl\r + l) 3 , 



64n 6 



0(n 



(A.4) 



r l/2( r _ 1 )3/2( r+1 )2^ r 2_ 1 ) 
64^ 



+ 0(n 



-8<\ 



and (X^pr)v = 5^/4 + 0(n- 8 ) for all other 1 < i, j < r. 

A straightforward computation that focuses on the lower-right two-by-two block of X^P^ 
shows that the matrix ^JX^P^ + | that appears in the expression for S n (3.3) has r — 1 
eigenvalues equal to 1 + 0(n~ 8 ) and one eigenvalue equal to 



(m 2 + 2)(2r + l) 2 «. 
1 + c n - c n ^ A + 0(n- 8 ) 



r 2 (r + l) 1 



(A.5) 



2n 2 v " 16n 4 

Similarly ■JXj n P^ h — | has r — 1 eigenvalues equal to 0(n -8 ) and one eigenvalue equal to 

(A.6) 



(m 2 + 2)(2r + l) 2 



2n 2 



Substituting these results into the expression for S n in (3.3) we obtain 

S n = c n (1 - logc n ) + c n logc n — + + - + 0(\ogn/n 8 ) . 



(A.7) 
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